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Abstract 

This paper presents the result of the investigation of chaotic osciUator synchroniza- 
tion. A new approach for detecting of synchronized behaviour of chaotic oscihators 
has been proposed. This approach is based on the analysis of different time scales 
in the time series generated by the coupled chaotic oscillators. This approach has 
been applied for the coupled Rossler and Lorenz systems. 
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1 Introduction 

Synchronization of chaotic oscillators is one of the fundamental phenomena 
of nonlinear dynamics. It takes place in many physical [1-5] and biological [6, 
7] processes. It seems to play an important role in the ability of biological 
oscillators, such as neurons, to act cooperatively [8-10]. 
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There are several different types of synclironization of coupled chaotic os- 
cillators which have been described theoretically and observed experimen- 
tally [11-14]. The complete synchronization (CS) implies coincidence of states 
of coupled oscillators Xi(i) = X2(t), the difference between state vectors of 
coupled systems converges to zero in the limit t ^ oo, [15-18]. It appears if 
interacting systems are identical. If the parameters of coupled chaotic oscil- 
lators slightly mismathch, the state vectors are close |xi(i) — X2(^)| w 0, but 
differ from each other. 

Another type of synchronized behavior of coupled chaotic oscillators wiht 
slightly mismatched parameters is the lag synchronization (LS), when shifted 
in time the state vectors coincide with each other, Xi(t + r) = X2(t). When 
the coupling between oscillator increases the time lag r decreases and the 
synchronization regime tends to be CS one [19-21]. 

The generalized synchronization (GS) [22-24], introduced for drive-responce 
systems, means that there is some functional relation between coupled chaotic 
oscillators, i.e. X2(t) = F[xi(t)]. Finally, it is necessary to mention the phase 
synchronization (PS) regime. To describe the phase synchronization the in- 
stantaneous phase (f){t) of a chaotic continuous time series is usually intro- 
duced [11-14,25,26]. The phase synchronization means the entrainment of 
phases of chaotic signals, whereas their amplitudes remain chaotic and uncor- 
related. 

All synchronization types mentioned above are concerned with each other (see 
for detail [1,20,22]), but the relationship between them is not completely clar- 
ified yet. For each type of synchronization there are its own ways to detect 
the synchronized behavior of coupled chaotic oscillators. The complete syn- 
chronization can be displayed by means of comparison of system state vectors 
xi(t) and X2(t), whereas the lag synchronization can be determined by means 
of the similarity function [19] 

^2. ) ^ (|x2(t + r)-xi(t)p) 



If the lag synchronization regime takes place the similarity function S{t) has 
its minimum a = min,- S{t) = for r corresponding to the time shift between 
the state vectors ^ . 

The case of the generalized synchronization is more intricate because the func- 
tional relation F[-] can be very complicated, but there are several methods to 



^ It is clear, that for the case of the complete synchronization S{t) reaches minimum 
value cr = for T = 0. 
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detect the synchronized behavior of coTiplcd chaotic oscillators, such as the 
auxiliary system approach [27] or the method of nearest neighbors [22,28]. 

Finally, the phase synchronization of two coupled chaotic oscillators occurs if 
the difference between the instantaneous phases (f){t) of chaotic signals Xi^2(^) 
is bounded by some constant 

1 01 (t) -02(^)1 < const. (2) 

It is possible to define a mean frequency 

n=]im^ = {m, (3) 

t— >oo t 



which should be the same for both coupled chaotic system, i.e., the phase 
locking leads to the frequency entrainmcnt. It is important to notice, to obtain 
correct results the mean frequency Cl of chaotic signal x(t) should coincide with 
the main frequency Qq = 27r/o of the Fourier spectrum (for detail, see [29]). 
Unfortunately, there is no general way to introduce the phase for chaotic time 
series. There arc several approaches which allow to define the phase for "good" 
systems with simple topology of chaotic attractor (so-called "phase coherent 
attractor"), the Fourier spectrum of which contains the single main frequency 
/o. The example of attractor and Fourier spectrum of such "good" system is 
shown in the figure 1. 

First of all, the instantaneous phase (pit) can be introduced as an angle in 
polar coordinates on the (x, |/)-plane [19,30] 

= arctan— , (4) 



but for that all trajectories of chaotic attractor projection on the {x, |/)-plane 
should revolve around some origin. Sometimes, a coordinate transformation 
can be used to obtain a proper projection [13,30]. Another way to define 
the phase (f)(t) of chaotic time series x{t) is the constructing of the analytical 
signal [11,25] 

at)^xit)+jx{t)^Ait)e^t>{t)^ (5) 

where the function x{t) is the Hilbcrt transform of x{t) 

x{t) = -P-V. I YZ~dT (6) 
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(where P.V. means that the integral is taken in the sense of the Cauchy princi- 
pal value). The instantaneous phase (f)(t) is defined from (5) and (6). Moreover, 
the Poincare secant surface can be used for the introducing of the instanta- 
neous phase of chaotic dynamical system [11,25] 

Ht) = 27r ^ ~^ + 27rn, tn<t< tn+i, (7) 

where t„ is the time of the nth crossing of the secant surface by the trajectory. 
Finally, the phase of chaotic time series can be introduced by means of the 
continuous wavelet transform [31], but the appropriate wavelet function and 
its parameters should be chosen [32]. 

All these approaches give correct results for "good" systems with well-defined 
phase, but fail for oscillators with non-revolving trajectories. Such chaotic 
oscillators are often called as "systems with ill-defined phase". The phase 
introducing based on the approaches mentioned above for the system with 
ill-defined phase leads usually to incorrect results [29]. Therefore, the phase 
synchronization of such systems can be usually detected by means of the 
indirect indications [30,33] and measurements [34]. 

In this paper we propose a new approach to detect the synchronization between 
two coupled chaotic oscillators. The main idea of this approach consists in the 
analysis of the system behavior on different time scales that allows us to 
consider different cases of synchronization from the universal positions. Using 
the continuous wavelet transform [35-38] we introduce into consideration the 
time scales s and associated with them instantaneous phases (j)s{t). As we'll 
show further, if two chaotic oscillators demonstrate any type of synchronized 
behavior, in the time series Xi(t) and X2(t) generated by these systems there 
are necessarily correlated time scales s for which the phase locking condition 

1 0.1 W -0.2^1 < const (8) 
is satisfied. 

The structure of this paper is the following. In section 2 we discuss the con- 
tinuous wavelet transform and the method of the time scales s and associated 
with them phases 4>s{t) definition. In the section 3 we consider the case of the 
phase synchronization of two mutually coupled Rossler systems. We demon- 
strate the application of our method and discuss the relationship between our 
and traditional approaches. The next section 4 deals with the synchronization 
of two mutually coupled Rossler systems with funnel attractors. In this case the 
traditional methods for the phase introducing fail and there is no possibility 
to detect the phase synchronization regime, respectively. The synchronization 
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between systems can be revealed here only by means of the indirect measure- 
ments (see for detail [34]). We demonstrate the efficiency of our method for 
such cases. In the section 5 we consider application of our method for the 
unidirectional coupled Rossler and Lorenz systems when the generalized syn- 
chronization is observed. In the section 6 we discuss the correlation between 
different types of the chaotic synchronization. The final conclusion is presented 
in the section 7. 



2 Continuous wavelet transform 



Let us consider continuous wavelet transform of chaotic time series x{t) 

W{s,to)^ J x{t)rs,toit)dt, (9) 



^s,to(t) = 1=J C—^) . (10) 



where ips,to{'t) is the wavelet-function related to the mother-wavelet ^o(^) as 

X- 

y/s \ S 

The time scale s corresponds to the width of the wavelet function ■il)s,tQ{t)-, 
and to is shift of wavelet along the time axis, the symbol "*" in (9) denotes 
complex conjugation. It should be noted that the time scale s is usually used 
instead of the frequency / of Fourier transformation and can be considered as 
the quantity inversed to it. 

The Morlet-wavelet [39] 

1 

i^o{^) = -77= exp(ifio^) exp — - (11) 



has been used as a mother-wavelet function. The choice of parameter value 
Qq = 271 provides the relation s = 1/f between the time scale s of wavelet 
transform and frequency / of Fourier transformation. 

The wavelet surface 

W{s,to) = \W{s,to)\e^^'^''^ (12) 



describes the system's dynamics on every time scale s at the moment of time 
to- The value of |Vr(s, to)| indicates the presence and intensity of the time 



5 



scale s mode in the time series x{t) at the moment of time to. It is possible to 
consider the quantities 



E{s,to) = \W{s,to)\'' 



(13) 



and 



■0, 



(14) 



which are instantaneous and integral energy distributions on time scales, re- 
spectively. 

At the same time, the phase 0s(t) = arg W{s,t) is naturally introduced for 
every time scale s. It means that it is possible to describe the behavior of 
each time scale s by means of its own phase (j)s{t). If two interacting chaotic 
oscillators are synchronized it means that in time series Xi(t) and X2(t) there 
are scales s correlated with each other. To detect such correlation one can 
examine the condition (8) which should be satisfied for synchronized time 
scales. 



3 Phase synchronization of two Rossler systems 

Let us start our consideration with two mutually coupled Rossler systems with 
slightly mismatched parameters [25, 26] 

yi,2 = <^i,2a:;i,2 + ayi,2, (15) 

Zl,2 = P + ^1,2(2^1,2 - C), 

where a = 0.165,p = 0.2, and c = 10. The parameters — ojQ:k/S. determine 
the parameter mistuning, e is the coupling parameter. In [19] it has been shown 
that for these control parameter values and coupling parameter e — 0.05 the 
phase synchronization is observed. 

For this case the phase of chaotic signal can be easily introduced in the one of 
the ways (4)-(7) mentioned above, because the phase coherent attractor with 
rather simple topological properties is realized in the system phase space. 
The attractor projection on the (x, i/)-plane resembles the smeared limit cycle 
where the phase point always rotates around the origin (Fig. l,a). The Fourier 
spectrum S{f) contains the main frequency peak /o ~ 0.163 (see Fig. 1,6) 
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Fig. 1. (a) Phase coherent attractor and (6) spectrum of the first Rossler sys- 
tem (15). Coupling parameter e between oscillators is equal to zero 

which coincides with the mean frequency / = Cl/2Tr, determinated from the 
instantaneous phase dynamics (3). Therefore, there are no comphcations 
to detect the phase synchronization regime in the two coupled Rossler sys- 
tems (15) by means of traditional approaches. 

When the coupling parameter e is equal to 0.05 the phase synchronization 
between chaotic oscillators is observed. The phase locking condition (2) is sat- 
isfied and the mean frequencies Qi_2 are entrained. So, the time scales Sq — 6 
of both chaotic systems corresponding to the mean frequencies fli^2 should 
be correlated with each other. Correspondingly, the phases (psi,2{t) associated 
with these time scales s should be locked and the condition (8) should be 
satisfied. The time scales which are the nearest to the time scale Sq should be 
also correlated, but the interval of correlated time scales depends upon the 
coupling strength. At the same time, should be time scales which remain un- 
correlated. These uncorrelated time scales cause the difference between chaotic 
oscillations of coupled systems. 

The figure 2 illustrates the behavior of different time scales for two coupled 
Rossler systems (15) with phase coherent attractors. It is clear, that the phase 
difference (j)si{t) — (ps2{t) for scales Sq = 6 is bounded and, therefore, time 
scales So = 6 corresponding to the main frequency of Fourier spectrum /o are 
synchronized. It is important to note that wavelet power spectra (i?i,2(s)) are 
close to each other (see Fig. 2, a) and time scales s characterized by the large 
value of energy (e.g., s=5) which are close to the main time scale sq = 6.0 
are correlated, too. There are also time scales which aren't synchronized, like 
s = 3.0, s = 4.0, etc. (see Fig. 2,b). 

So, when two mutually coupled chaotic oscillators with phase coherent attrac- 
tors are considered, the traditional method based on the instantaneous phase 
0(i) of chaotic signal introducing and our approach lead to the equivalent 
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Fig. 2. (a) Wavelet power spectrum {E{s)) for the first (solid line) and the sec- 
ond (dashed line) Rossler systems (15). (6) The dependence of phase difference 
(t>si{t) — 4>s2it) on time t for different time scales s. The couphng parameter between 
oscillators is £ = 0.05. The phase synchronization for two coupled chaotic oscillators 
is observed 

results. 



4 Synchronization of two Rossler systems with funnel attractors 



Let us consider more complicated example when it is impossible to correctly 
introduce the instantaneous phase 4>{t) of chaotic signal x(t). It is clear, that 
for such cases the traditional methods of the phase synchronization detecting 
fail and it is necessary to use another techniques, e.g., like indirect measure- 
ments [34]. On the contrary, our approach gives correct results and allows to 
detect the synchronization between chaotic oscillators easily as before. 

To illustrate it we consider two non-identical coupled Rossler systems with 
funnel attractors (Fig. 3): 

il,2 = -<^l,2yi,2 - ^1,2 + £ix2,l - 3^1,2), 

yi,2 = t^i,2a:;i,2 + ayi,2 + £(^2,1 - ^1,2), (16) 

Zl,2 = P + ^1,2(2^1,2 - C), 



where £ is a coupling parameter, uji = 0.98, UJ2 = 1-03. The control parameter 
values have been selected by analogy with [34] as a = 0.22, p — 0.1, c = 8.5. 

In [34] it has been shown by means of the indirect measurements that for the 
coupling parameter value e = 0.05 the synchronization of two mutually cou- 
pled Rossler systems (16) takes place. Our approach based on the analysis of 
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the dynamics of different time scales s gives analogous results. So, the behavior 
of the phase difference (f)si{t) — (f)s2{t) for this case has been presented in the 
figure 4,6. One can see that the phase locking takes place for the time scales 
s — 5.25 which are characterized by the largest energy value in the wavelet 
power spectra {E{s)) (see Fig. 4, a). Thus, we can say that the time scales 
s = 5.25 of two oscillators are synchronized with each other. It is important to 
note that the other time scales (e.g., s = 4.5,6.0 et. al.) remain imcorrelated. 
For such time scales the phase locking has not been observed (see Fig. 4,6). 

It is clear, that the mechanism of the synchronization of coupled chaotic os- 
cillators is the same in both cases considered in the sections 3 and 4. The 
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Fig. 3. (a) Phase picture and (6) power spectrum of the first Rossler system (16) 
oscillations. Coupling parameter e is equal to zero 
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Fig. 4. (a) The normalized energy distribution in wavelet spectrum {E(s)) for the 
first (the solid line denoted as "1") and the second (the dashed hne denoted as "2") 
Rossler systems (16); (6) the phase difference (psiit) — (t)s2{t) for two coupled Rossler 
systems. The value of coupling parameter has been selected as e = 0.05. The time 
scales s = 5.25 are correlated with each other and the synchronization has been 
observed 
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synchronization phenomenon is caused by the existence of time scales s in 
system dynamics correlated with each other. Therefore, there is no reason to 
divide considered synchronization examples into different types. 



5 Generalized synchronization of Rossler and Lorenz systems 

Let us consider another type of synchronized behavior, so-called the gener- 
alized synchronization. As oscillator samples the Rossler (drive) and Lorenz 
(response) systems have been selected. The equations of the drive system 
xi = {xi,yi,zi)'^ are 

xi = -ujyi - zi, 

yi^ujxi + ayi, (17) 
ii ^p + zi{xi - c), 

and the response system X2 = (2:2, 2/2, 2^2)"^ is given by 
X2 = -a{x2 - 1/2), 

y2 = ru(t) - 2/2 - u{t)z2, (18) 
Z2 = u{t)y2 - bz2. 

The drive signal u{t) and control parameter values have been selected by 
analogy with [23] as u{t) — Xi + yi + zi, p — 2, c = 4, u = 1, a — 1, a — 10, 
r = 28, b = 2.666. It has been analytically shown (see [23]) that the response 
system (18) is asymptotical stable for arbitrary drive signals u{t) and arbitrary 
initial conditions, and, therefore the generalized synchronization always occurs 
although drive and response systems are completely different. Obviously, the 
generalized synchronization regime can be also detected by means of other 
numerical or experimental methods (e.g., the axillary system approach). 

The chaotic attractors and corresponding them power spectra S{f) of unidi- 
rectional coupled Rossler and Lorenz systems are shown in the figures 5 and 
6, respectively. One can see in the Fourier spectrum of the Lorenz system the 

presence of the peaks corresponding to the frequencies of the Rossler system 
oscillations (see Fig. 5,b and 6,6). Therefore, the time scales s of coupled sys- 
tems can be correlated with each other if chaotic oscillators are synchronized 
on these time scales. 
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Fig. 6. (a) Chaotic attractor and (b) spectrum of the Lorenz (response) system (18). 



The dependencies of phase difference (f)si{t) — (f>s2{t) on time for different time 
scales s are shown in the figure 7. It is clear, that there is the range of scales 
(approximately s = 5 -r- 7) the behavior of which is synchronized. At the same 
time there are time scales which are not synchronized, as in the case of the 
phase synchronization (see. e.g.. Fig. 2). It is important to note, despite the 
fact that there is the frequency /r ~ 0.33 in the spectra of both Rossler 
and Lorenz systems (see Fig. 5,5 and 6,6), the time scales sr — — 
3 of coupled oscillators are not synchronized (Fig. 7,6). So, the presence of 
equivalent frequencies / in the power spectra of interacting systems doesn't 
mean the obligatory synchronization on the time scale Sf = 1// corresponding 
to this frequency. 

Thus, the generalized synchronization of two unidirectionally coupled abso- 
lutely different system is revealed as the synchronous dynamics of several time 
scales on which the phases (j)s{t) are locked. One can see that different types of 
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Fig. 7. (a) The normalized energy distribution in wavelet spectrum {E(s)) for the 
Rossler (17) (the solid line denoted as "1") and Lorenz (18) (the dashed line de- 
noted as "2") systems. Each energy distribution has been normalized on its own 
coefficient; (6) the phase difference (f)si{t) — (t)s2{t) for coupled systems, when the 
generalized synchronization regime takes place. The time scales s = 5, s = 6, s = 7 
are correlated with each other whereas the other time scales (e.g., s = 3, s = 8, 
etc.) aren't synchronized 

the synchronization such as PS and GS are quite similar to each other when 
the time scale behavior of interacting systems is analysed. So, the time scales 
of Rossler and Lorenz systems in the generalized synchronization regime be- 
have themselves equivalently to the cases of the phase synchronization of two 
Rossler systems (15), although for the considered generalized synchronization 
case we can not correctly introduce the instantaneous phase (f)(t) for neither 
Rossler nor Lorenz systems (see Fig. 5 and 6) at all. Obviously, it is neces- 
sary to consider the correlation between different types of synchronization and 
transitions from one of them to another one. This topic will be discussed in 
the next section. 



6 From unsynchronized behavior to complete synchronization regime 

It has been shown in [19] that there is certain relationship between PS, LS 
and CS for chaotic oscillators with slightly mismatched parameters. With the 
increase of coupling strength the systems undergo the transition from unsyn- 
chronized chaotic oscillations to the phase synchronization. With a further 
increase of coupling the lag synchronization is observed. The next increasing 
of the coupling parameter leads to the decreasing of the time lag and both 
systems tend to have the complete synchronization regime. 

Let us consider the dynamics of different time scales s of two nonidentical 
mutually coupled Rossler systems (16). If there is no phase synchronization 
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between the oscillators, then their dynamics remain uncorrelated for all time 
scales s. The figure 8 illustrates the dynamics of two coupled Rossler systems 
when the coupling parameter e is enough small {e = 0.025). The power spec- 
tra {E(s)) of wavelet transform for Rossler systems differ from each other 
(Fig. 8, a), but the maximum values of the energy correspond approximately 
to the same time scale s in both systems. It is clear, that the phase difference 
4>siit) — (f)s2it) is not bounded for all time scales (see Fig. 8,b). It means that 
there are no time scales s correlated with each other and two Rossler systems 
aren't synchronized at all. 
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Fig. 8. (a) The normalized energy distribution in wavelet spectrum {E(s)) for the 
first (the solid line denoted as "1") and the second (the dashed line denoted as "2") 
Rossler systems; (b) the phase difference (j)si{t) — 4>s2{t) for two coupled Rossler 
systems. The value of coupling parameter has been selected as £ = 0.025. There is 
no phase synchronization between systems 

As soon as any of the time scales of the first chaotic oscillator becomes 
correlated with the other one of the second oscillator (e.g., when the cou- 
pling parameter increases), the phase synchronization occurs (see Fig. 4 in 
the section 4). It is clear, that the time scales s characterized by the largest 
value of energy in wavelet spectrum {E{s)) become correlated first. The other 
time scales remain uncorrelated as before. The phase synchronization between 
chaotic oscillators leads to the phase locking (8) on the correlated time scales 
s. 

When the parameter of coupling between chaotic oscillators increases, more 
and more time scales become correlated and one can say that the degree of 
the synchronization grows. So, with the further increasing of the coupling 
parameter value (e.g., e — 0.07) in the coupled Rossler systems (16) the time 
scales which were uncorrelated before become synchronized (see Fig. 9,6). One 
can see that the time scales s = 4.5 arc synchronized in comparison with the 
previous case {e = 0.05, Fig. 4,6) when these time scales were uncorrelated. 
The number of time scales s demonstrating the phase locking increases, but 
there are nonsynchronized time scales as before (e.g., the time scales s — 3.0 
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and s — 6.0 remain still nonsynchronized) . 
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Fig. 9. (a) The normalized energy distribution in wavelet spectrum {E{s)) for the 
first (the solid line denoted as "1") and the second (the dashed line denoted as "2") 
Rossler systems; (6) the phase difference (t)si{t) — (t)s2{t) for two coupled Rossler 
systems. The value of coupling parameter has been selected as e = 0.07. 



The arising of the lag synchronization [19] between oscillators means that all 
time scales are correlated. Indeed, from the condition of the lag-synchronization 
xi{t — r) ~ X2{t) one can obtain that Wi{s,t — r) ~ W2{t,s) and, therefore, 
(psi{t — r) ~ (f)s2{t). It is clear, in this case the phase locking condition (2) is 
satisfied for all time scales s. E.g., when the coupling parameter of chaotic os- 
cillators (16) becomes enough large (s = 0.25) the lag synchronization of two 
coupled oscillators appears. In this case the power spectra of wavelet transform 
coincide with each other (see Fig. 10, a) and the phase locking takes place for 
all time scale s (Fig. 10,6). It is important to note that the phase difference 
0si(O — 4>s2{t) is not equal to zero for the case of the lag synchronization. It 
is clear that this difference depends on the time lag r. 

The next increasing of the coupling parameter leads to the decreasing of the 
time lag r [19]. Both systems tend to have the complete synchronization regime 
Xi{t) ~ X2{t), therefore the phase difference (f)si{t) — (t>s2{t) tends to be a zero 
for all time scales. 

So, it is clear, that PS, LS and CS are naturally interrelated with each other 
and the synchronization type depends on the number of synchronized time 
scales. At the same time, the relationship between PS and GS is not clear 
in detail. There are several works [1, 20] dealing with the problem, how GS 
and PS arc correlated with each other. E.g., in [20] it has been reported that 
two unidirectional coupled Rossler systems can demonstrate the generalized 
synchronization while the phase synchronization has not been observed. This 
case allows to be easily explained by means of the time scale analysis. The 



14 




a b 



Fig. 10. (a) The normalized energy distribution in wavelet spectrum {E(s)) for 
the Rossler system; (b) the phase difference (psiit) — <ps2{t) for two coupled Rossler 
systems. The value of coupling parameter has been selected as £ = 0.25. The lag 
synchronization has been observed, all time scales are synchronized 

equations ol Rossler system are 

xi = -(jJiyi - zi, 
iji = ujixi + ayi, 
Zi =p+ Zi{Xi -c) 

X2^ -UJ2y2- Z2 + e{xi-X2), (19) 

y2 = i^2X2 + ay2, 

Z2 ^P + Z2{X2 - C), 



where xi = {xi,yi, Zi)'^ and X2 = {x2, 1/2, ^2)'^ are the state vectors of the first 
(drive) and the second (response) Rossler systems, respectively. The control 
parameter values have been chosen as uji = 0.8, LJ2 = 1-0, a = 0.15, p = 0.2, 
c = 10 and e — 0.2. The generahzed synchronization takes place in this case 
(see [20] for detail). Why it is impossible to detect the phase synchronization in 
the system (19) although the generalized synchronization is observed becomes 
clear from the time scale analysis. 

Let us consider Fourier spectra of coupled chaotic oscillators (see Fig. 11). 
There are two main spectral components with frequencies /i = 0.125 and 
/2 = 0.154 in these spectra. The analysis of behavior of time scales has shown 
that both the time scales Si = I//1 = 8 of coupled oscillators corresponding 
to the frequency /i and time scales close to Si are synchronized while the 
time scales S2 = I//2 — 6.5 and close to them don't demonstrate synchronous 
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Fig. 11. Fourier spectra for(a) the first (drive) and (6) the second (response) Rosier 
systems (19). The couphng parameter is e = 0.2. The generahzed synchronization 
takes place 

behavior (Fig. 12,b). 
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Fig. 12. (a) The normalized energy distribution in wavelet spectrum {E{s)) for the 
first (the solid line denoted as "1") and the second (the dashed line denoted as 
"2") Rossler systems. The time scales pointed by means of arrows correspond to 
the frequencies /i = 0.125 and /2 = 0.154, respectively; (6) the phase difference 
(t>si{t) — (i^s-iit) for two coupled Rossler systems. The generalized synchronization 
has been observed 

The source of such behavior of time scales become clear from the wavelet 
power spectra {E{s)) of both systems (see Fig. 12, a). The time scale Si of 
the drive Rossler system is characterized by the large value of energy while 
the part of energy falling on this scale of the response system is quite small. 
Therefore, the drive system dictates its own dynamics on the time scale Si to 
the response system. The contrary situation takes place for the time scales S2 
(see Fig. 12, a). The drive system can not dictate its dynamics to the response 
system because the part of energy falling on this time scale is small in the first 
Rossler system and enough large in the second one. So, time scales S2 are not 
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synchronized. 

Thus, the generahzed synchronization of the unidirectional coupled Rossler 
systems appears as the time scale synchronized dynamics, as before another 
synchronization types do. It is also clear, why the phase synchronization hasn't 
been observed in this case. One can see from the Fig. 11 that the instantaneous 
phases 0i,2(^) of chaotic signals Xi.2(t) introduced by means of traditional 
approaches (4)-(7) are determined by the both frequencies /i and /2, but 
only the spectral components with the frequency /i arc synchronized. So, the 
observation of instantaneous phases 0i,2(^) doesn't allow to detect the phase 
synchronization in this case although the synchronization of time scales takes 
place. 

Thus, one can see that there is a close relationship between different types 
of the chaotic oscillator synchronization. Unfortunately, it is not clear, how 
one can distinguish the phase synchronization^ and the generalized synchro- 
nization using only the results obtained from the analysis of the time scale 
dynamics and what kind of the relationship between these synchronization 
types is. We suppose that GS and PS are often practically equivalent (when, 
of course, it is possible to define correctly the instantaneous phase of chaotic 
signal by means of traditional technique). Nevertheless, this problem should 
be separately investigated later. 



7 Conclusion 

Summarizing this work we would like to note several principal aspects. Firstly, 
the traditional approach for the detecting of the phase synchronization based 
on the introducing of the instantaneous phase (f){t) of chaotic signal is suit- 
able and correct for such time series which are characterized by the Fourier 
spectrum with the single main frequency /q. In this case the phase 0^0 as- 
sociated with the time scale sq corresponding to the main frequency /o of 
the Fourier spectrum coincides approximately with the instantaneous phase 
(f){t) of chaotic signal introduced by means of the traditional approaches (see 
also [32]). Indeed, as the other frequencies (the other time scales) don't play 
a significant part in the Fourier spectrum, the phase 0(t) of chaotic signal is 
close to the phase 0so(^) of the main spectral frequency /o (and the main time 
scale So, respectively). It is obvious, that in this case the mean frequencies 
/ = (0(t))/27r and fso — (0so(i))/27r should coincide with each other and with 



^ We mean here the phase synchronization between chaotic oscillators takes place 
if the instantaneous phase (p{t) of chaotic signal may be correctly introduced by 
means of (4)-(7) and the phase locking condition (2) is satisfied. 
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the main frequency /o of the Fourier spectrum (see also [29]) 

f = Lo = k. (20) 

If the chaotic time scries is characterized by the Fourier spectrum without the 
main single frequency (like the spectrum shown in the Fig. 3,6) the traditional 
approaches (4)-(7) fail. It is clear that one has to consider the dynamics of 
the system on all time scales, but it is impossible to do by means of the 
instantaneous phase 0(t) introduced as (4)-(7). On the contrary, our approach 
based on the continuous wavelet transform can be used for both types of 
chaotic signals. 

Secondly, our approach can be easily applied to the experimental data because 
it doesn't require any a-priori information about the considered dynamical 
systems. Moreover, in several cases the influence of the noise can be reduced 
by means of the wavelet transform (for detail, see [38,40,41]). We believe that 
our approach will be useful and effective for the analysis of physical, biological, 
physiological etc. data, such as [8,31,32]. 

Finally, it is important to note that analysis of the system dynamics on the dif- 
ferent time scales based on the continuous wavelet transform allows to consider 
the different types of behavior of coupled oscillators (such as the complete syn- 
chronization, the lag synchronization, the phase synchronization, the general- 
ized synchronization and the nonsynchronized oscillations) from the universal 
position. It is clear, that the number of synchronized time scales determines 
the type of behavior uniquely. Probably, the quantitative characteristic of the 
synchronization measure can be introduced. This method (with insignificant 
modifications) can also be applied to dynamical systems synchronized by the 
external (e.g., harmonic) signal. 



8 Acknowledgments 

We express our appreciation to George A. Okrokvertskhov, Alexander V. 
Kraskov and Professors Vadim S. Anishchenko and Tatyana E. Vadivasova 
for valuable discussions. We thank also Svetlana V. Eremina for the support. 



18 



References 



[I] Parlitz U., Junge L., Lauterborn W., Experimental observation of phase 
synchronization, Phys. Rev. E 54 (2) (1996) 2115-2117. 

[2] Tang D.Y., Dykstra R., Hamilton M.W., Heckenberg N.R., Experimental 
evidence of frequenct entrainment between coupled chaotic oscillators, Phys. 
Rev. E 57 (3) (1998) 3649-3651. 

[3] AUaria E., Arccchi F.T., Garbo A.D., Mcucci R., Synchronization of homoclinic 
chaos, Phys. Rev. Lett. 86 (5) (2001) 791-794. 

[4] Ticos CM., Rosa E., Pardo W.B., Walkcnstein J.A., Monti M., Experimental 
real-time phase synchronization of a paced chaotic plasma discharge, Phys. Rev. 
Lett. 85 (14) (2000) 2929. 

[5] Rosa E., Pardo W.B., Ticos CM., Wakenstein J.A.,Monti M., Phase 
synchronization of chaos in a plasma discharge tube. Int. J. Bifurcation and 
Chaos 10 (11) (2000) 2551-2563. 

[6] Tass P. A. at al., Detection of n:m phase locking from noisy data: Application 
to magnetoencephalography, Phys. Rev. Lett. 81 (15) (1998) 3291-3294. 

[7] Anishchenko V.S., Balanov A.G., Janson N.B., Igoshcva N.B., Bordyugov G.V., 
Entrainment between heart rate and weak nonlinear forcing, Int. J. Bifurcation 
and Chaos 10 (10) (2000) 2339-2348. 

[8] Elson R.C at al.. Synchronous behavior of two coupled biological neurons, Phys. 
Rev. Lett. 81 (25) (1998) 5692-5695. 

[9] Rulkov N.F., Modeling of spiking-bursting neural behavior using two- 
dimensional map, Phys. Rev. E 65 (2002) 041922. 

[10] Tass P. A. et al.. Synchronization tomography: A method for three-dimensional 
localization of phase synchronized neuronal populations in the human brain 
using magnetoencephalography, Phys. Rev. Lett. 90 (8) (2003) 088101. 

[II] Pikovsky A., Rosenblum M., Kurths J., Phase synchronisation in regular and 
chaotic systems. Int. J. Bifurcation and Chaos 10 (10) (2000) 2291-2305. 

[12] Anishchenko V.S., Vadivasova T.E., Synchronization of self-oscillations 
and noise-induced oscillations, Journal of Communications Technology and 
Electronics 47 (2) (2002) 117-148. 

[13] Pikovsky A., Rosenblum M., Kurths J., Synhronization: a universal concept in 
nonlinear sciences, Cambridge University Press, 2001. 

[14] Anshchenko V.S., Astakhov V., Neiman A., Vadivasova T., Schimansky- 
Geier L., Nonlinear Dynamics of Chaotic and Stochastic Systems. Tutorial and 
Modern Developments, Springer- Verlag, Heidelberg, 2001. 

[15] Pecora L.M., Carroll T.L., Synchronisation in chaotic systems, Phys. Rev. Lett. 
64 (8) (1990) 821-824. 



19 



[16] Pecora L.M., Carroll T.L., Driving systems with chaotic signals, Phys. Rev. A 
44 (4) (1991) 2374-2383. 

[17] Murali K., Lakshmanan M., Drive-response scenario of chaos syncronization in 
identical nonlinear systems, Phys. Rev. E 49 (6) (1994) 4882-4885. 

[18] Murali K., Lakshmanan M., Transmission of signals by synchronization in a 
chaotic van der polduffing oscillator, Phys. Rev. E 48 (3) (1994) R1624-R1626. 

[19] Rosenblum M.G., Pikovsky A.S., Kurths J., From phase to lag synchronization 
in coupled chaotic oscillators, Phys. Rev. Lett. 78 (22) (1997) 4193-4196. 

[20] Zheng Z., Hu G., Generalized synchronization versus phase synchronization, 
Phys. Rev. E 62 (6) (2000) 7882-7885. 

[21] Taherion S., Lai Y.C., Observability of lag synchronization of coupled chaotic 
oscillators. Phys. Rev. E 59 (6) (1999) R6247-R6250. 

[22] Rulkov N.F., Sushchik M.M., Tsimring L.S., Abarbanel H.D.I., Generalized 
synchronization of chaos in directionally coupled chaotic systems, Phys. Rev. E 
51 (2) (1995) 980-994. 

[23] Kocarev L., Parlitz U., Generalized synchronization, predictability, and 
equivalence of unidirectionally coupled dynamical systems, Phys. Rev. Lett. 
76 (11) (1996) 1816-1819. 

[24] Pyragas K., Weak and strong synchronization of chaos, Phys. Rev. E 54 (5) 
(1996) R4508-R4511. 

[25] Rosenblum M.G., Pikovsky A.S., Kurths J., Phase synchronization of chaotic 
oscillators, Phys. Rev. Lett. 76 (11) (1996) 1804-1807. 

[26] Osipov G.V., Pikovsky A.S., Rosenblum M.G., Kurth J., Phase synchronization 
effect in a lattice of nonidentical rossler oscillators, Phys. Rev. E 55 (3) (1997) 
2353-2361. 

[27] Abarbanel H.D.I., Rulkov N.F., Sushchik M., Generalized synchronization of 
chaos: The auxiliary system approach, Phys. Rev. E 53 (5) (1996) 4528-4535. 

[28] Pecora L.M., Carroll T.L., Heagy J.F., Statistics for mathematical properties of 
maps between time series embeddings, Phys. Rev. E 52 (4) (1995) 3420-3439. 

[29] Anishchenko V.S., Vadivasova T.E., Journal of Communications Technology 
and Electronics 49 (1). 

[30] Pikovsky A., Rosenblum M., Osipov G., Kurths J., Phase synchronization of 
chaotic oscillators by external driving, Physica D 104 (4) (1997) 219-238. 

[31] Lachaux J. P. at al., Studying single-trials of the phase synchronization activity 
in the brain, Int. J. Bifurcation and Chaos 10 (10) (2000) 2429-2439. 

[32] Quiroga R.Q., Kraskov A., Krcuz T., Grassberger P., Perfomance of different 
synchronization measures in real data: a case study on electroencephalographic 
signals, Phys. Rev. E 65 (2002) 041903. 



20 



[33] Pikovsky A.S., Rosenblum M.G., Kurths J., Synchronization in a population of 
globally coupled chaotic oscillaors, Europhysics Letters 34 (3) (1996) 165-170. 

[34] Rosenblum M.G., Pikovsky A.S., Kurths J., Locking-based frequency 
measurement and synchronization of chaotic oscillators with complex dynamics, 
Phys. Rev. Lett. 89 (26) (2002) 264102. 

[35] Daubechies I., Ten lectures on wavelets, SIAM, 1992. 

[36] Kaiser G., A friendly guide to wavelets, Springer Verlag, 1994. 

[37] Torresani B., Continuous wavelet transform, Paris: Savoire, 1995. 

[38] Koronovskii A. A., Hramov A.E., Continuous wavelet analysis and its 
applications, Moscow, Fizmatlit, 2003. 

[39] Grossman A. and Morlct J., Decomposition of hardy function into square 
integrable wavelets of constant shape, SIAM J. Math. Anal. 15 (4) (1984) 273. 

[40] Torrence C, Compo G.P., A practical guide to wavelet analysis. Bulletin of the 
American Meteorological 79 (1) (1998) 61-78. 

[41] Gusev V.A., Koronovskiy A. A., Hramov A.E., Application of adaptive wavelet 
bases to analysis of nonlinear systems with chaotic dynamics, Technical Physics 
Letters 29 (18) (2003) 61-69. 



21 



